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Quantum communication devices, such as quantum repeaters, quantum memories, or quantum 
channels, are unavoidably exposed to imperfections. However, the presence of imperfections can 
be tolerated, as long as we can verify such devices retain their quantum advantages. Benchmarks 
based on witnessing entanglement have proven useful for verifying the true quantum nature of 
these devices. The next challenge is to characterize how strongly a device is within the quantum 
domain. We present a method, based on entanglement measures and rigorous state truncation, 
which allows us to characterize the degree of quantumness of optical devices. This method serves 
as a quantitative extension to a large class of previously-known quantum benchmarks, requiring no 
additional information beyond what is already used for the non-quantitative benchmarks. 



I. INTRODUCTION 

Quantum communication is the transmission of infor- 
mation via quantum states. Such states can exhibit cor- 
relations which are far stronger than anything allowed 
by classical theory. These quantum correlations enable 
us to perform novel communication tasks, such as the 
teleportation of quantum states [1] or the distribution 
of secret cryptographic keys through public channels [2]. 
To support the transfer and storage of quantum infor- 
mation, our communication infrastructure must incorpo- 
rate new quantum devices such as quantum repeaters, 
quantum teleportation systems, quantum memories, and, 
most generally, quantum channels. Unfortunately, un- 
avoidable imperfections in these devices can degrade the 
embedded quantum information. Although we can deal 
with some degradation by employing classical postpro- 
cessing, if the imperfections become too strong we can 
lose the quantum advantage altogether. Thus, it is im- 
portant to distinguish whether our communication de- 
vices are truly functioning in a quantum way. 

Given a quantum state as input, a classical device is 
one which uses a measure and prepare strategy (see Fig. 
la) . This strategy consists of a measurement on the state, 
storage or transmission of the measurement result classi- 
cally, and preparation of a new output state based on this 
classical information. Any device which can be shown to 
be inconsistent with this strategy is said to operate in the 
quantum domain. From this basic idea, we can develop 
benchmarks which tell us whether or not a device is in 
the quantum domain. For example, one approach is to 
use the fidelity as a figure of merit. If the fidelity between 
the output state and the input state, averaged over some 
ensemble of test states, is higher than the best fidelity 
achievable by a measure and prepare strategy, then the 
device is in the quantum domain. Unfortunately, this 
approach has some practical drawbacks. To calculate 
the fidelity accurately, we need to either perform (exper- 
imentally costly) full tomography on the output states or 
make unverified assumptions about them (e.g., that they 
are Gaussian) . As well, the optimal classical fidelity must 



be found anew for each type of test ensemble, and it has 
only been found for a few cases [3-9]. Moreover, many of 
these fidelity benchmarks require testing with an infinite 
ensemble of states, which is not possible in experiment. 

Another approach to the quantum benchmarking prob- 
lem relies on witnessing entanglement, since measure and 
prepare channels are mathematically equivalent to the so- 
called entanglement breaking channels [10]. If we verify 
that a channel or device preserves entanglement, when 
acting on a subsystem, then we also know that the de- 
vice operates in the quantum domain. In fact, this can 
even be done with a minimum of experimental resources. 
For one, this approach does not actually require the use 
of entangled states; instead, we can use 'effective' entan- 
glement [11-14]. As well, only a few test states and a few 
different measurements suffice to witness entanglement. 
This experimental simplicity, without the need for un- 
verified assumptions, makes entanglement-based bench- 
marks very practical. This approach is especially useful 
for testing optical devices with continuous degrees of free- 
dom [15, 16], where the fidelity-based techniques cannot 
be easily used. Entanglement-based benchmarks have al- 
ready been designed which use two [14, 17] or more [18] 
coherent test states, two squeezed test states, and even 
mixed test states [19]. 

These benchmarks (fidelity- or entanglement-based) do 
not give the full story so far. Currently, they only ad- 
dress the simple question: quantum or classical? There 
is a much richer realm of possibilities within the quan- 
tum domain. A device could be in the quantum domain, 
yet operate very inefficiently, consuming more resources 
than another device which performs the same task. To 
capture this difference in performance, we need quanti- 
tative extensions of existing benchmarks. Such measures 
will enable us to characterize and compare different de- 
vices within the quantum domain. For entanglement- 
based benchmarks, a natural extension is to consider how 
well devices preserve entanglement, by making use of en- 
tanglement measures. Such an approach was recently 
studied in [20] , calling this notion the quantum through- 
put. This quantitative extension required no additional 
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The remainder of this paper is organized as follows. 
In Sec. II, we outline the main notions of entanglement- 
based benchmarking, both qualitative and quantitative. 
Then, in Sec. Ill, we present our new method for cal- 
culating quantitative benchmarks. We apply our scheme 
to a variety of scenarios and report on the results in Sec. 
IV. We conclude the paper with Sec. V, while detailing 
some proofs in Appendices A and B. 



II. BENCHMARKING QUANTUM OPTICAL 
DEVICES USING ENTANGLEMENT 



A. Qualitative benchmarking 



FIG. 1: (Color online) a) Device implementing a measure and 
prepare strategy: A measurement is made on the input state, 
the result is stored or transmitted classically, and a new quan- 
tum state is prepared based on this measurement result, 
b) Example distributions of test states: Entanglement-based 
benchmarks have been designed for ensembles consisting of 
two or more coherent states on a ring, squeezed and an- 
tisqueezed vacuum, as well as for generalizations to mixed 
states. 



experimental resources beyond the (already limited) re- 
quirements of previous entanglement-based benchmarks. 

The results of [20] are an important first step, but there 
remains much room for improvement. For one, the 'quan- 
titative domain' should be as faithful as possible with the 
quantum domain. In other words, if a device preserves 
the entanglement of some test state, then we should aim 
to give a non-trivial estimate of how much entanglement 
remains. We may not be able to estimate the amount 
exactly, but it is certainly non-zero. Unfortunately, the 
methods from [20] only give non-trivial estimates for a 
small fraction of the quantum domain (see Sec. II), lead- 
ing to overly stringent device requirements. Clearly there 
is a need for more work on the theoretical side. Addition- 
ally, [20] considers only one type of benchmarking situ- 
ation, involving two coherent test states. It is desirable 
to have quantitative extensions for many of the other 
entanglement-based benchmark schemes, and we need to 
develop new tools to handle these cases. 

In this paper, we present a new method to quanti- 
tatively extend existing entanglement-based benchmarks 
for continuous variable optical devices. This new method 
is very successful, giving much stronger results than were 
previously achievable, often faithful over nearly all of the 
quantum domain. The method relies on projecting opti- 
cal states to large finite subspaces and rigorously bound- 
ing the resulting truncation errors using homodyne mea- 
surement results. The strength of our approach can also 
be improved simply by devoting more computational re- 
sources to the task. The method is quite general and can 
be applied to many of the existing entanglement-based 
benchmarking situations, requiring no additional experi- 
mental resources. 



Here, we briefly review some basic notions of 
entanglement-based benchmarking schemes. The inter- 
ested reader is referred to the references for further de- 
tails. These schemes are based around the fact that the 
set of measure and prepare channels is exactly the same 
as the set of entanglement breaking channels [10]. If we 
represent a device by a superoperator A, then it is entan- 
glement breaking if the state [1 <g> A] (p) is separable for all 
p. To disclude this behaviour, we must find a state which 
remains entangled after interaction with the device. To 
this end, consider the following entangled state: 
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where Alice's subsystem is d-dimensional (using the com- 
putational basis) and Bob's subsystem is an optical 
mode. The \tpk) come from some fixed ensemble of 
nonorthgonal test states. The nonorthogonality of these 
states is encoded in Alice's reduced density matrix: 



(m\ p A \n) = - (V> m | rj} n ) ■ 



(2) 



Alice's subsystem is kept isolated, while Bob's subsystem 
is sent through the device. This leads to the final state 
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Since Alice's subsystem was kept isolated, the reduced 
density matrix of this final state is still given by Eq. (2). 

Now, Alice measures her subsystem using the projec- 
tive measurement {| JZq. When Alice measures the 
state \j), Bob's subsystem is projected, in general, into 
an infinite-dimensional mixed state 
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We call these 'conditional states,' since they are condi- 
tioned on Alice's measurement outcome. Bob measures 
the conditional states using a balanced homodyne detec- 
tion setup. This allows him to determine the expectation 
values (x), (ft), and variances 



Var(x) = (x 2 ) - (x) 2 
Var(p)=(p 2 )-(p} 2 



(5) 
(6) 
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FIG. 2: (Color online) a) Entanglement-based benchmarking: 
Is there any entanglement remaining after part of the state 
|^ ent ) passes through the device? How much remains? 
b) 'Effective' entanglement: Any source which produces a 
classical mixture of the test states \iph) can be pictured as 
internally preparing the entangled state in Eq. (1) and per- 
forming a projective measurement on subsystem A. 
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FIG. 3: (Color online) Example benchmark results using two 
coherent states (T = 0.448). The upper line was found using 
the EVM method; any data below this line is in the quantum 
domain. The lower line represents the limits of the quantita- 
tive domain using the methods of [20]. There is clearly room 
for improvement on the quantitative benchmarks. 



of the canonical quadratures operators, 

x = —j=(a) + a), p = — ^(a 1 — a) 
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This procedure is repeated many times to get useful mea- 
surement statistics for each of the conditional states. The 
goal is for Alice and Bob to verify the continued presence 
of entanglement in the final state pab- The given infor- 
mation (homodyne measurement results and knowledge 
of pa) turns out to be sufficient for this task. 

In fact, if we use a so-called 'prepare and measure' 
scheme 1 [11, 14], it is not necessary to actually use an 
entangled state to test if a device is entanglement break- 
ing. Instead, we can use a classical mixture of the test 
states \ipk): an d ask whether the device would have pre- 
served entanglement if this mixture represented one sub- 
system of the entangled state Eq. (1). Any device which 
produces such a classical mixture can be thought of this 
way (see Fig. 2). Such entanglement is called effective 
entanglement, and this way of thinking offers a great ex- 
perimental advantage. 

To map out the quantum domain, we need to parama- 
terize typical measurement results. For the quadratures, 
the action of any device on a test state can be described 
using the induced loss and noise. If the initial test states 
\ipu) have means (x in ) k , (pi n ) k and variances Var fc (£ m ), 
Varfe(p in ), then the loss is parametrized by 1 — T, where 
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Not to be confused with measure and prepare strategies. 



and the excess noise is introduced through 

Var k (x out ) = Var k (x in ) + V ex (x). (9) 

Here, we assume that all conditional states exhibit the 
same loss and excess noise, and also that these quanti- 
ties are symmetric with respect to x and p. This is done 
merely for simplicity of presentation, and the entangle- 
ment verification or quantification methods do not rely 
on this choice. 

One tool which has proved quite useful for the task 
of entanglement verification is the so-called expectation 
value matrix (EVM) [14], which has been successfully 
used in all of the benchmarking scenarios mentioned 
above [17-19] (with suitable substitutions in Eq. (1)). 
The EVM method allows us to determine, for a given 
set of test states, which measurement results provide 
evidence that a device operates the quantum domain. 
An example result of the EVM method is shown in 
Fig. 3, for the case of testing with two coherent states 
|^)e{|a>,|-a». 



B. Quantitative benchmarking 

The main focus of this paper is in quantitative exten- 
sions to the above device testing scenario. The overall 
goal is to quantify how well an optical device preserves 
quantum correlations (the quantum throughput) by es- 
timating the amount of entanglement, with respect to 
some measure £ , that remains in the output state pab 
of Eq. (3). We will develop the entanglement quantifi- 
cation method in general, but later use the negativity 
[21-23] in concrete examples. The available information 
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is the same as what is used for the qualitative benchmark- 
ing, namely homodyne measurements on the conditional 
states f? B and knowledge of pa . Since this information is 
tomographically incomplete, the best we can do is to give 
a lower bound on the entanglement, by minimizing over 
all states (Jab which have the same expectation values 
[24]: 
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min £{ctab) 
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= Tr(p B M k ) 
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The measurement operators may be something di- 
rectly measured in experiment (such as the quadratures) 
or something more abstract (such as projections onto 
pure states). 

First results for this quantification problem were pre- 
sented in [20] . The underlying challenge is that the states 
in question are, in general, infinite-dimensional, making 
quantification inherently difficult. The basic strategy in- 
troduced to tackle this was to consider a projection of 
the infinite-dimensional mode subsystem down to a more 
manageable finite-dimensional subspace. This is more 
than just a practical simplification. It has been shown 
that, for any state pab, the presence of entanglement can 
be verified by searching within some finite subspace [25] . 
In fact, for any entanglement measure with the strong 
monotonicity property, the entanglement content of the 
projected state gives a lower bound to the entanglement 
of the true state p AB . 

To see this, let P be a projector onto the desired mode 
subspace, and denote the (unnormalized) projected state 
by ' 

p P := (l A ®P)p AB (l A ®P). (11) 
Then we have (with Q := t B — P and p = Tr(pp)) 
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For practical purposes, it is convenient to absorb the pref- 
actor p and work directly with the unnormalized states 
(0 < Tr pp < 1). As such, we restrict to entanglement 
measures which are well-defined for unnormalized states 
and where, for any p £ [0, 1] and r > 0, 

P £{t) > £{pr). (13) 

In this case, we end up with the final inequality 

£{pab) > £{pp)- (14) 



The above constraints must be kept in mind when choos- 
ing an entanglement measure. 

To go along with this projection, we must also relax 
the measurement constraints in Eq. (10). Importantly, 
even if we know an expectation value of pab exactly, we 
may not be able to determine the corresponding expec- 
tation value of pp exactly. However, we should try to 
constrain the expectation values of the projected state 
to be as close as possible to those of the state. If we 
understand the system well enough, we can bound how 
much the expectation values can vary, ideally restricting 
them to some convex set C which depends on the known 
measurement results. If we can find such constraints, 
then we can perform a minimization over all projected 
states ap compatible with the known constraints about 
the actual pp (and hence about pab)' 
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(15) 



Finally, by solving the optimization of Eq. (15) , we will 
arrive at a rigorous lower bound on the entanglement of 
the infinite-dimensional output state pab, 



£(Pab) > £min(pp)- 



(16) 



To summarize, instead of finding the minimal entangle- 
ment over infinite-dimensional states with exactly known 
expectation values, we find the minimal entanglement 
over finite states with inexact expectation values. This 
latter minimization gives a lower bound to the former. 
Any choice of projection requires us to find appropriate 
constraints on the expectation values, and the tighter the 
constraints we can find, the stronger the results will be. 

But what subspace should be projected onto? In [20], 
which considered only the testing scenario involving two 
coherent states (\tpi) £ {\a) , |— a)}), the choice is made 
to project onto a two-qubit subspace, the smallest sub- 
space which can show entanglement. In this case, the 
operators M k are projections onto the 'largest' eigenvec- 
tors of the final mode states p l B (i — 0,1). Constraints 
on these eigenvectors can be determined from the homo- 
dyne measurements using some bounds from the litera- 
ture [26, 27]. A typical result of this approach is pre- 
sented in Fig. 3. The quantum domain, found using the 
EVM method, is shown in the same figure. 

Although the two qubit projection works, it is clear 
that there are some drawbacks, since the quantitative 
domain does not come close to matching the quantum 
domain. First, by projecting into the smallest possible 
subspace, we may be cutting out important information. 
Intuitively, we expect most information about our states 
to be confined to a small subspace, but a two-qubit sub- 
space appears to be too small. Second, the employed 
bounds are only tight for very low noise values; this weak- 
ens the minimization at higher noise levels. Finally, the 
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bounds are applicable only to a specific two-qubit sub- 
space, and not readily generalized to higher dimensions. 

In the following section, we propose a modified version 
of this quantification scheme which proves to be much 
more successful. We still make use of the projection to a 
finite subspace, but we choose a different type of subspace 
which can have significantly larger dimension. For this, 
we will derive new, rigorous bounds which allow us to 
make much tighter approximations on our state. This 
new method does not require any additional information, 
and unlike the previous quantification scheme, is readily 
generalized to other benchmarking scenarios. 



The coefficients Ck are positive, but otherwise unknown. 
We have 



Tr(rn) = ^ kr k = X! Ck = n > 



and 



fc=0 k=l 



N 



This leads to 



(21) 



Tr(r N n) = ^c k = n N . (22) 
fe=l 



III. PROJECTING ONTO FOCK STATES 
A. Rigorous bounds on the projection errors 

Perhaps the most natural choice of projection in the 
infinite-dimensional Hilbert space of an optical mode is a 
projection onto the energy eigenstates, namely the Fock 
states {|n)}„_ . Intuitively, the larger the cutoff N is, 
the closer the projected state pp should be to the full 
state pab, and for N — > oo, the two states are the same. 
Although we will use the finite projection trick of [20] 
to find entanglement bounds, the projection onto Fock 
states requires a much different set of tools for its analysis 
than the previously considered two-qubit projection. 

Most importantly, we make the observation that it is 
not necessary to send N to infinity. As long as the cut- 
off dimension N is large compared to the mean photon 
number n, a state and its projected form will be nearly 
identical. Why is this so? Because high Fock levels con- 
tribute more to the state's energy than low Fock levels. 
If we know n, we can put strong bounds on the contri- 
butions from high Fock levels. We make this intuition 
concrete with the following lemma. 

Fock cutoff lemma. Let r be the state of a single optical 
mode, with mean photon number Tr(rn) = n < oo. Let 
tjv be the (unnormalized) projection of r up to the Fock 
level TV > 0, 



N 

TN = J2(k\r\k)\k)(k\ 

k=0 



(17) 



and denote the expectation value of t/v with respect to 
n by 



Tr(Tjvn) =: n N . 
Then the following inequality holds: 



Tr(r) - Tr(rjv) < 



n — n N 
N+ 1 ' 



(18) 



(19) 



Proof: For k > 0, parametrize the diagonal entries of r 
by 



Tfe = (k\ T \k) 



Tr(r)-Tr(r N )= £ | 

k=N+l 



< 



ELfe 
N + l 
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N 

n — tin 
'N+l 



(23) 
(24) 

(25) 
(26) 



Finally, note that this upper bound can actually be sat- 
urated, for example by the state 



-=^1-^10X01 + — |A + 1)(A+1|. (27) 

□ 

This lemma provides an upper bound on how different 
the trace of a state is from the trace of its projected 
form. Since the mean photon number can be inferred 
from homodyne measurements, 



Tr(rn) = Tr r 



- P 



(20) 



(28) 



we can easily apply this lemma to the present bench- 
marking scenario. In this case, the quantity tin will be a 
free parameter which is a linear combination of diagonal 
elements from t/v- It is only constrained by the obvious 
bounds < njv < n. Although it would be valid to re- 
place the right hand side of Eq. (19) with the strictly 
numeric quantity n/{N + 1), it is more useful to use the 
bound as is, since it links the constraints on Tt(tat) and 
un- Most importantly, since n is fixed, the bound can be 
made arbitrarily tight by increasing the cutoff N. Along 
the same lines, we can derive similar constraints (also 
depending on n and N) for all the other important ex- 
pectation values. Together, these constraints will form 
the basis of our new entanglement quantification method, 
outlined in the next subsection. 

For instance, if a state r has some known values of (x) 
and (p), then the expectation values of the same oper- 
ators with respect to the truncated state tjv cannot be 
entirely arbitrary. In fact, if N is large enough, the trun- 
cated state's quadratures have to be quite close to the 
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infinite-dimensional state's quadratures. We can make 
this argument explicit with the following bound, involv- 
ing the ladder operator a = ^ (x + ip) : 



with 



\Tr( T a)-Tr(T N a)\<e, (29) 



Vfa-njvXTtM-Tty-jv-i)), (30) 



where tjv_i is the truncation of the state r up to Fock 
level N—l. This constraint can be straightforwardly de- 
rived by expanding the annihilation operator in the Fock 
basis and applying the triangle and Cauchy-Schwarz in- 
equalities, as well as appealing to the positivity of r. The 
full derivation can be found in Appendix A. Although 
the quantity Tr(ra) is not directly measureable, it can 
be found from the expectation values of the quadratures 
x and p with respect to r. As well, the bound e is not 
fixed, since it depends on other parameters. However, by 
the Fock cutoff lemma, it has a maximum size 



e < 



n 
N 



(31) 



Thus, it is clear that the bound (29) can also be made 
tighter by increasing N. In practice, we will use the 
stricter adaptive constraint, Eq. (30). 

The other important operators for our problem are x 2 
and p 2 . Instead of directly working with these, it is useful 
to consider their difference, 



V 



;t2 



(32) 



Although this operator is itself an observable, the re- 
quired expectation values can be readily found from the 
available homodyne data. The original operators can be 
recovered using the relations 



P 2 = 



d + 2h + l 
-d + 2n + 



(33) 
(34) 



In a similar manner as before, we can bound the differ- 
ence between (^dtj for the measured infinite-dimensional 
state and for its truncation: 



with 
5:=2yf( 



Tr(rd) - Tr(T N d) 
a — n N ) [{n 



(35) 



n N -2) + (Tr(r) - Tr(r A r_ 2 ))]. 

(36) 

The proof of this bound uses similar arguments as the e 
bound (see Appendix B for details). 

The final source of information in our benchmarking 
scheme is the reduced density matrix pa , which enters in 
expectation values of the form 



To see how to include this type of information, let tab 
be a bipartite state, with reduced density matrix ta = 
I^bCtab)- As above, let P be a projector on the B sub- 
system up to Fock level N, and let Q = 1b — P (which 
implies QP = PQ = 0). We also let 



:= (1 a ®P)tab(Ia®P), 



(38) 



and similarly define tq using Q. Using the given rela- 
tions, we find 

T A =Tt b {tab(1a ®[P + Q] 2 )) (39) 
=Ty b (tab(1a <8 Pf) + Ty b (tab(1a ® Q) 2 )- (40) 

But 

Tvb(tab{1a®P) 2 ) = Tr B ({l A ®P)T AB (lA®P)), (41) 

(and similarly for the Q term), so that 

T A = Tr B (T N )+Tr B (T Q ). (42) 

Since tn and tq are positive matrices, so are their re- 
duced forms t n a ■= Ttb(tjv) and tqa := Tr b (tq). 
Dropping tqa, we are left with 



TA > TNA- 



(43) 



Tr(cTAB |i)01 ® : 



(j I Pa \i) 



(37) 



Although it may not be obvious at this point, Eq. (43) 
can be used to tightly constrain expectation values like 
Eq. (37) for the projected state. As with the other 
bounds, increasing the cutoff N will strengthen the con- 
straints. We will elaborate further about how to use this 
inequality below. The constraints found in Eqs. (29-30), 
(35-36), and (43) form the basis for determining a rigor- 
ous lower bound to a state's entanglement content. We 
explore how to do this in practice in the next subsection. 



B. Finding the minimal entanglement with 
semidefinite programming 

These new constraints can be directly applied to the 
problem of quantitative benchmarking. We know the ex- 
pectation values of l,a, and d for the conditional states 
p l B of pab- These serve as the operators M fc found in 
the optimizations (10) and (15). As well, we have the 
entries of the reduced density matrix pa- Using the 
above constraints, we can restrict these same expecta- 
tion values for the projected state pjq (we will use the 
subscript N instead of P from now on to emphasize the 
Fock state projection). We search through the projected 
subspace for the state with the minimal entanglement. 
The above bounds constrain more strictly where we can 
search, forming the convex regions C(p B ,Mk) and C(pa) 
in the finite-dimensional optimization (15). The larger 
we make the cutoff N, the tighter the constraints will 
be, and asymptotically, these bounds become equality 
constraints, making the optimizations (10) and (15) the 
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same. Finally, we also enforce the obvious positivity 
((Tat > 0) and trace (Tr(o-jv) < 1) constraints. 

To solve the finite-dimensional optimization (15), we 
turn to semidefinite programming. In principle, the ma- 
trices representing the operators x and p in the Fock ba- 
sis are infinite-dimensional. This can lead to problems 
when working with them numerically. However, for a 
state which only has support up to Fock level N, there 
is no cause for concern. Such expectation values can be 
implemented by using truncated versions of the quadra- 
ture operators themselves. Indeed, if xn is the (bounded) 
projection of x up to Fock level N, then 



Taking the determinant, we get 



Tt(t n x) = Tt(t n x n ). 



(44) 



Expectation values of the quadratures are simply finite 
linear combinations of elements from r^r. The same idea 
holds as well for higher powers of the quadratures, as in 
the operator d. 

Conveniently, the above constraints can all be written 
in the form of matrix inequalities. For instance, for any 
r which is positive, the constraint in Eqs. (29-30) may 
be recast as 



n — tin Tr(ra) — Tr(rjva) 

Tr(rat)-Tr(rjvat) Tr(r) - Tr(rjv-i) 



> 0. (45) 



Similarly, we can rewrite the second order constraint, 
Eqs. (35-36), as the constraint 



4(n - n N ) 
Tr(rd) — Tr(rjvrf) n 



Tr(rd) - Tr(r N d) 
n N -2 + Tr(r) - Tt(t N - 2 ) 



> 0. 

The original constraints follow from considering the pos- 
itivity of the determinants. These matrix inequality con- 
straints can be applied for each conditional state p B . 
Finally, we have the constraint 



PA > Pna- 



(47) 



To see how this applies in a benchmarking situation, we 
note that the output state pab can be written in the 
block form 



PAB 



Poo Poi P02 

PlO Pll Pl2 
P20 P21 P22 



(48) 



where the blocks pa correspond to the conditional states 
p l B . The (i, j) element of pa is exactly equal to Tr(pij)/d. 
We can similarly decompose p^ into a block form where 
the matrices are replaced by their projected forms pf^ 
up to Fock level N. The element of pNA is analo- 
gously equal to Tr(p^)/d. The inequality (47) allows us 
to constrain the trace of each block pfj. To see this, con- 
sider a scenario where subsystem A is two dimensional. 
Then matrix inequality (47) becomes 



Tr(poo) 
Tr (/9io) 



Tr(p^ ) 



Tr(poi) 
Tr(pii) 



Tr(pK) 



> 0. (49) 



|Tr(p i)-Tr(^)| < 7 , 



(50) 



where 



7 := V (Tr(poo) 



Tr(^ ))(Tr( Pll ) 



Tr(pfO)- (51) 



This constraint is therefore very similar to the other ones, 
and it should be clear that we can make it stronger by 
increasing the cutoff N. We should point out that if sub- 
system A has more than two dimensions, the inequality 
(47) becomes even stronger since it also forces higher or- 
der determinants to be positive. 

We now have the full set of constraints on p^ for the 
minimization, but we have yet to fix our objective func- 
tion, i.e. the actual entanglement measure. As men- 
tioned earlier, in order for the finite projection to work 
as desired, the measure should have the strong mono- 
tonicity property and also allow us to absorb prefactors 
and work with unnormalized states. As well, the measure 
(or a lower bound on it) should be computable within the 
framework of semidefinite programming. To this end, we 
choose the negativity [21-23], though this is not the only 
possible choice. Although the negativity may not have as 
nice an operational interpretation as some other entan- 
glement measures, its main strength is that it is efficiently 
computable. For a Hermitian matrix, the negativity can 
be written as the following semidefinite program 2 : 



(52) 



A/"(t) = min Tr r_ 

s.t. t Ta = t + — t_ , 

T+ > 0, T_ > 0, 

where t Ta is the partial transpose of r. 



IV. RESULTS 

The new quantification method can be used to extend 
many existing entanglement-based benchmarks. We will 
present results for the following ensembles of test states: 
a) two coherent states [17, 26], b) two squeezed states 
[19], and c) three coherent states [18]. All optimizations 
were done in Matlab with the solver SDPT3 [28] and the 
frontend YALMIP [29], using a Dell Optiplex 760 desktop 
PC with a 3 GHz dual-core processor and 4 GB of RAM. 

First, our results for the two coherent state benchmark 
are shown in Fig. 4, where we have used a cutoff at 
N = 20 and fixed the loss parameter as T = 1.0. We 
can see that the results are very faithful, being non-zero 



2 For comparison purposes, it is important to note that that in 
[20] the negativity is rescaled by a factor of 2 from the formula 
used here. 
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FIG. 4: (Color online) Entanglement bounds for the two co- 
herent state benchmark, using our new quantification scheme 
(with N = 20), assuming no loss (T = 1.0). The dashed 
black line represents the boundary of the quantum domain, 
found using the EVM method. Variances lower than this line 
correspond to devices within the quantum domain, whereas 
variances higher than this line are classical [17] (note that the 
axes are inverted in order to see the results better). Direct 
comparison with the results of [20] (Fig. 2a) shows a dramatic 
improvement in the entanglement bounds at higher variances. 



Classical 
domain 




Var Q (x) 



FIG. 5: (Color online) Entanglement quantification for 
squeezed/antisqueezed state benchmarks, using initial squeez- 
ing magnitude r = 0.35 and cutoff N = 20 (cf. Fig. 1 of [19]). 
The bottom left region is forbidden by the uncertainty rela- 
tion. The quantum domain lies between the forbidden region 
and the upper dotted line (found with EVM method). The 
quantitative domain, shown using level curves of the minimal 
negativity, is very faithful in this parameter regime; the zero 
contour and the quantum domain boundary coincide within 
the numerical resolution. For much larger values of Varo (x) , 
the two domains begin to diverge due to reduced numeri- 
cal accuracy. The small 'sawtooth' gap near the unphysical 
boundary is due to the finite resolution of sampled points. 



throughout nearly all of quantum domain. For other val- 
ues of T, where the quantum domain is smaller, the re- 
sults are similarly faithful. Only in the extreme regions 
of low overlap/low variance and high overlap/high vari- 
ance are there small differences between the quantitative 
domain and the quantum domain. Both of these cases 
correspond to high values of n, where the constraints are 
the weakest. But we expect very little effective entangle- 
ment can be found in these regions anyway. The region 
where (a\— a) — > 1 has very little initial entanglement 
(since the states are nearly identical), and there will be 
even less after subsystem B passes through the device. 
On the other hand, in the region where (a\— a) — > 0, 
Pa becomes diagonal and our source replacement trick is 
useless. Although there may be a lot of real entangle- 
ment, we cannot verify any effective entanglement with a 
diagonal pa- In principle, the quantitative domain could 
be extended by raising the cutoff N. Since in the limit 
N — > oo, the optimization problem (15) converges to 
(10), we conjecture that any mismatch between the quan- 
titative domain and the quantum domain can be resolved 
by using sufficiently large computational resources. 

Another important class of entanglement-based bench- 
marks uses squeezed and antisqueezed vacua as test 



states. To simplify the display of results, we restrict our- 
selves to devices whose operation is phase independent, 
so that the changes in the long/short quadratures are the 
same for both test states, regardless of their orientation 
[19]. In this case, the relevant parameters for the quan- 
tum domain are simply the initial squeezing magnitude 
r and the output variances of one test state, Varo (x) and 
Varo(p). Example results for this benchmark type are 
shown in Fig. 5, for the parameters r = 0.35 and N = 20. 
We see that the quantitative domain is very faithful with 
the quantum domain in the given parameter region. For 
much larger values of the variance than shown in Fig. 5, 
the entanglement bounds eventually decay. Again, this is 
due to the chosen cutoff N = 20 not being large enough 
to counter higher noise values. 

Finally, the last example of entanglement-based bench- 
marking we will consider is for three coherent states dis- 
tributed symmetrically on a ring. Benchmarks with mul- 
tiple coherent states are known to lead to larger quantum 
domains than those with just two coherent states [18]. 
The magnitude \a\ of the coherent states is used as an 
input parameter, while loss and excess noise paramatrize 
the output. The results of applying our quantification 
procedure to this scenario are shown in Fig. 6, using 
N = 15 (this is lower than in the previous cases because 
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FIG. 6: (Color online) Entanglement quantification for the 
three coherent state benchmark, with no loss (T = 1.0) and 
cutoff TV = 15 (cf. Fig. 3 of [18]). Variance values below 
the dashed black line are in the quantum domain (note again 
that the axes are flipped to aid visualization of the results). 
Although the quantitative domain is larger than the case with 
two coherent states, it does not cover the entire quantum 
domain due to the smaller cutoff and higher energies involved. 

the extra dimension on subsystem A also needs to be ac- 
comodated within our available computation resources). 
We find that the quantitative domain is larger than for 
two coherent states. However, there is a small but notica- 
ble gap between the quantitative domain and the known 
quantum domain. 

To explore the strength of our scheme, we consider the 
following hybrid of the optimizations (10) and (15): 

{min S((Tn) 
<TJV>0 
s.t. Tr(cr pf ® M fc ) = Tr(p 4 B M fe ) . 
Tr(a N \i)(j\®l B ) = (j\p A \i) 

(53) 

where, as above, the operators Mk are taken from the 
set {l,a, d}. In this optimization problem, we search 
through states in the projected subspace (for some fixed 
value of N), but we force the expectation values to be 
exactly those of the infinite-dimensional state pab ■ The 
result of this optimization will necessarily be an upper 
bound on the minimal entanglement of pab- In Fig. 7, 
we compare the upper bound found by this hybrid opti- 
mization (using N = 15) with the results shown in Fig. 
6 by considering a section at amplitude |a| = 0.2. We 
see that our result is fairly tight to the upper bound for 
most variance values. We also see that there is very little 
negativity to be found (less than 0.03) for Var(x) > 1. 
The mismatch between the quantitative domain and the 



0.25 




0.5 0.6 0.7 0.8 0.9 1 1.1 1 .2 1 .3 1 .4 1 .5 

Var(x)=Var(p) 

FIG. 7: (Color online) Upper and lower bounding the minimal 
entanglement. The lower curve is a slice through the surface in 
Fig. 6 at \a\ = 0.2 and the upper curve corresponds to results 
of the hybrid optimization (53). We see that the minimal 
entanglement must be quite small after Var(x) as 1. Higher 
values of the cutoff N would be needed to have small enough 
error bounds to detect this entanglement. 



quantum domain again seems to be due to the combina- 
tion of low entanglement and weaker constraints at high 
noise values. If desired, the quantitative domain could 
likely be extended by employing more computational re- 
sources. 



V. CONCLUSION 

We have outlined a new method for quantitatively ex- 
tending entanglement-based benchmarks. The method 
provides rigorous lower bounds on the entanglement con- 
tent of test states after partial action by an optical device. 
The method, which requires no additional information, 
leads to very good agreement with many of the known 
quantum domains, providing a huge improvement over 
the quantification scheme of [20] (though we should note 
that the results presented here do not supersede those of 
[20] , which appears to be a better choice for devices with 
very low levels of imperfection). Most importantly, the 
method presented here enables us, for the first time, to 
quantitatively characterize the quantumness of realistic 
optical devices. 
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Appendix A: Derivation of first order constraint 



Appendix B: Derivation of second order constraint 



Here, we detail how to arrive at the bounds given in 
Eqs. (29-30). Let r be a bounded positive-scmidcfinitc 
operator and let tn be the projection of r onto the 
subspace given by the first N + 1 Fock states, tm = 
J2k=o (M T \ty' We seek to bound the magnitude of the 
difference 



Tr(Ta) — Tr(rAra) = (k\ (r — tn)cl \k) 

k=0 

oo 

= ^y/k{k\{T-T N )\k-l). (Al) 



fc=0 



For simplicity, denote general matrix elements by 
(m|r|n) =: r m n , and diagonal entries by T m m —: r m . 
Then, from the triangle inequality, 



|Tr(ra) - Tt{t n ci)\ < ^ ^\ r k,k-i\ 

k=N+l 

From the positivity of r, we must have 

I'm 7~n ^ | 7~m . n \ 

Therefore, 

oo 

|Tr(ra) - Tr(T N a)\ < ^ ^\f^k\J T k-\ 

k=N+l 



(A3) 



(A4) 



Using the Cauchy-Schwarz inequality on this sum, we 
find 



|Tr(Ta)-Tr(T W a)| < 



\ 



kTk )[ H rfc -! ■ 

\k=N+l / \k=N+l / 

(A5) 

The first bracketed term is the difference in mean pho- 
ton number between t and t/v, while the second is the 
difference in trace between r and its truncation after the 
level TV — 1. Finally, we arrive at the desired bound: 



|Tr(ra) - Tr(r N a)\ < \/{n- ti n )(Tt:(t) - Tr(rjv-i))- 

(A6) 



In this appendix, we outline how to arrive at the bound 
on the difference operator found in Eqs. (35-36). We 
begin with 



Tr(ra t2 ) - Tr(r w at2) = (k\ (r - r N )o) 2 \k) 

k=0 

oo 



J2 Vk + Wk + 2r k , k +2- 
(Bl) 



k=N-l 



Similarly, 

oo 

Tr^a^-Tr^a 2 ) = ^ Vk + 1 y/k + 2r* k+2 . (B2) 



k=N-l 



(A2) Therefore, with d = a) 2 + a 2 , 



Tr(rd) - Tr(TjvO 



2 J2 Vfc + iVk + 2Re[T fc , fe+2 ] 

fe=JV-l 
oo 

<2 Vk + lVk + 2\T k , k+2 \ . 



k=N-l 



(B3) 



As in Appendix A, we can use the positivity of r and the 
Cauchy-Schwarz inequality to get 



Tr(rd) - Tr(r N d) 



< 2 



\ 



£ (k + 2)Tt +2 ) Y. (* + !)'>]• (B4) 



\k=N-l 



\k=N-l 



Finally, by rewriting the right hand side as 

2a/ (n - n N ) [(n - n N -2) + (Tr(r) - Tr(r w _ 2 ))] = 5, 

(B5) 

we arrive at the desired bound. 
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